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1 Introduction 


In this paper we propose a new method of calculation of generating functions of the generalized 
Chebyshev polynomials associated with root systems of simple Lie algebras. The method is 
illustrated by calculations of the generating functions of the two-variables Chebyshev polyno¬ 
mials associated with the simple Lie algebras A 2 , C 2 and G 2 - 

The Chebyshev polynomials in several variables are a natural generalization of the classical 
polynomials in one variable |T]. They are used in different areas of mathematics (for example, 
in discrete analysis, in approximation theory j2], in linear algebra 0 . 0 . in representation 
theory m-m) and in physics [8] - [12]. 

The classical Chebyshev polynomials of the first kind Tn{x) is defined by the formula 

Tn{x) = Tn{cos(f)) = cosncf), n = 0,1,... (1) 

They satisfy the following three-term recurrence relation 

Tn+l{x) = 2xTn(x) - Tn-l{x), U > 1, (2) 

with the initial conditions 

To(a;) = 1, Ti(x)=x. (3) 

The simplest d-dimensional generalization of the Chebyshev polynomials is reduced to the 
product of one-dimensional polynomials Tn^Xi) 

d 

rn(x) = ]^r„,(xi), n = (ni, ...,nd), x = (xi,..., x^). 

However, much more interesting multidimensional analogues of the polynomials ([T]) which are 
associated with roots systems of simple Lie algebras na - [I9]. Recall their definition following 
mainly the work [16| . 

Let d-dimensional Euclidean space is divided into regularly spaced cells representing 
convex polyhedra. Such a division is similar to splitting the axis at intervals of periodicity 
of the function cos0. We assume that there are many hyperplanes that divide any of these 
polyhedra on a finite number of small congruent polyhedra, similar to the original one (in 
the classical one-dimensional case each hyperplane degenerates into the point that divides the 
interval of periodicity in two equal parts). Such polyhedra are called foldable. The reflection 
of each polyhedron with respect to any hyperplane translates it into another polyhedron of 
splitting. Such a cellular structure of folding polyhedra are defined by a corresponding group 
of reflections. 

Let L be a simple Lie algebra and R be a reducible system of roots. A system of roots is a 
set of vectors in d-dimensional Euclidean space supplied with a scalar product (.,.) which 
is completely determined by a basis of simple roots a*, i = 1,.., d and by a group of reflections 
of R called Weyl the group W{R). Generating elements of the Weyl group Wi, i = I, ..,d act 
on any vector x ^ E‘^ according to the formula 


2(x, ttj) 

WiX = X — - - -tti- 

Specifically, if x = Oj we obtain from (0]) tCjCTj = —a*, 
corresponds the coroot 

V 

a = - --. 

(a, a) 


( 4 ) 


To any root a from the system R 
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For the basis of the simple coroots i = 1,c? one can dehne the dual basis of fundamental 
weights Aj, z = 1,d 

(Ai, aj) = 

(here and in what follows the dual space E'^* is identihed with E^). The bases of roots and 
weights are related by the linear transformation 


a* = 



a, 




(5) 


where C is the Cartan matrix of the Lie algebra L. 

The set of weights A satisfying the inequality (A, ai) > 0, where z = 1, d, are the simple 
roots, is called the dominant Weyl chamber Cq. The images of Cq obtained by action on it by 
all elements of the Weyl group are called the Weyl chambers or fundamental domains. Thus 
the Weyl group acts on the set of Weyl chambers transitively. Weyl chambers play the role of 
convex polyhedra in the partitioning of the space E^ which was described above. A reflection 
of a Weyl chamber with respect to any hyperplane reduces to composition of the action of 
elements of the Weyl group and translation group. The semidirect prod^t of a Weyl group W 
and a comnmtative translation group T is called an affine Weyl group W of a root system R. 
The group W is just a symmetry group of the cell structure of the space E'^. 

A function which can be considered as a possible candidate for ad- dimensional general¬ 
ization of cos n(j), i.e., multivariable Chebyshev polynomial, must be invariant under the action 
of the Weyl group and have the property of periodicity. Following [TB] - |19| . we dehne the 
periodical function in d variables (orbit function in terms of the work jTB]) by the formula 


$„(0) = ^ 
UlgW 


( 6 ) 


which is obviously W - invariant ^u,n{(p) = ‘^*n(0), Vzn G W. Here n is expressed in the basis 
of fundamental weights {Aj} and 4> is expressed in the dual basis of coroots {a/} 

d d 

n = ^ UiXi riiE Z, ^ (j)i e [0,1). 

i=l i=l 


The function $11(0) dehned by the relation (jS]), with non-negative integer rii from the n = 
(zzi,..., rz^), up to a normalization gives the d-variable Chebyshev polynomials of fhe 

hrst kind. 

For the function $ 11 ( 0 ) we have the following multiplication rule 


*hk*hs — 'y ^ *h.y,k;+s- 

w£W 


(7) 


This formula allows us to obtain the recurrent relations for the corresponding polynomials of 
several variables. 

The function $ 11 ( 0 ) from (jS]) can be represented in the form 


$n(0) 


\W\ d 

\ 27ri fki(^) hii 

- 

k=l 
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where fki{(p) are some linear functions of (0i,..,0d) and \W\ is the order of the Weyl group 
(that is the number of its elements). Let us introduce the following diagonal matrix 


M 


diag 


27ri fli(0) Di 

e i=i 


27ri Y, f|W|i(0) Hi 

e 


( 8 ) 


In view of (jS]) we have 


<hn(0) = tr M. 


The matrix M can be rewritten as a product of diagonal matrices 


d 

M = 

k=l 


( 9 ) 


( 10 ) 


where 


Mk = diag g27rif|w|k(0)) = g27nAk^ ^ f\w\k{(l>)) ■ (H) 

Now let us associate to each Mk the following matrix 

Rk = {I\w\ — PkMk) (12) 

where I\w\ is the unit |hL| x |hL|-matrix and pk is a real parameter. This matrix Rk is different 
from the resolvent matrix for Mk in the location of the spectral parameter. Since all the matrices 
Rk are diagonal, they commute with each other. Obviously, the following relation is hold 


1 d°>^ 

^ (i?,j 


Ukl dpj)'" 


(13) 


Pk=0 


It follows from ([9]), (ITOl) and (lT3l) that 

4>„(<A) = = tr(Mr'..AO = 


Qni+..+na 


(tr(i?p,..i?pj) 


ni\..nd\ dpi^-.dpn'^ 

The following function of the parameters {pi} 

^L;Pd = '^dRpi-Rpd = '^nr,..,nA4})pT-p7, 

ni..nd>0 


Pl=--=Pd=0 


(14) 


(15) 


can be considered as a (non-normalized) generating function of Chebyshev polynomials in d 
variables, although is still defined in terms of exponent . The expression of 

tr(i?pj..i?p^) in terms of the variables Xt defined as 


2—1 d—i 

Xj = $0^(0), Bi = (C^, liC^), (16) 

is not too cumbersome. This will be illustrated in the following examples of the Lie algebras 
A 2 , C *2 and G 2 - 

Calculation of the generating function of Chebyshev polynomials in several variables of the 
second kind requires a slight modification of the method described above. 


4 







Classical Chebyshev polynomials of the second kind Un{x) is defined by the relation 


Unix) 


sin (n + 1)0 
sin 0 


n > 0. 


( 17 ) 


They satisfy the same recurrence relations ([ 2 ]) as the polynomials of the first kind with different 
initial conditions 

Uoix) = 1, Uiix) = 2x. (18) 

Generalization of polynomials flTT)) for the case of several variables based on the Weyl character 
formula |21] . |22) . from which it follows 


t/n(0) 


^ detw 
togW 

X) det w ’ 

togW 


(19) 


where dettc = (—1)^1’") and iiw) is the minimal number of Weyl group generating elements 
Wi required for expressing to as a product of Wi. p is the Weyl vector which is equal to the 
half-sum of the positive roots. Using the formula ([5]) the Weyl vector can be expressed in the 
basis of fundamental weights. 

The function 

^ detw 

which stands in the numerator of flT^ . can be represented as the difference between the expres¬ 
sions and 


^ n 




E 

det w=l 


g27ri(ii;(n+p),</)) 


E 

iogW, det w=—l 


g27ri(w(n+p),0) 


Further calculation repeats the scheme given above ([H]) - flT^ for each of the functions 
separately. In the considered case we must use the relation 




d“ 


Till..rial d°ipi..d“dpj. 


(tr(K 


'pi'-^Pd ^P1 




( 20 ) 


Pi = --=Pd=0 


instead of <!>„ given by (fTTl) . In the formula (|20|) the matrices R^. are used to represent functions 
More details are presented in the examples below. 

To obtain Chebyshev polynomials of the second kind we must in accordance with the formula 

d 

(IT^ divide the function by the singular element <F“, where 0 = (0, ..,0). Thus, the 
generating function of Chebyshev polynomials of the second kind has the form 


„ tr(fl;...fi* - fi-..fl-) 

PI--,Pd <1,^ 


( 21 ) 


As in the case of the first kind polinomials, transition to the variables Xi = Uei{4>) in Fp^ 
does not require cumbersome calculations. 

In the second part of this work we calculate the generating functions of the two variables 
Chebyshev polynomials associated with the systems of simple roots for the Lie algebras A 2 , C 2 
and G 2 . 
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2 The case of Lie algebra A 2 

2.1 The generating function for the polynomials of the first kind 

In this section we obtain the well known two variables Chebyshev polynomials associated with 
the root system of the Lie algebra A 2 by the technique presented in the Section 1. These 
polynomials were introduced by Koornwinder in [13], and their generating functions were found 
in the work [2^ by another method. 

The reduced root system R of the Lie algebra A 2 includes two fundamental roots ai, a 2 , 
positive root ai + 02 and the reflections of these roots. The action of the Weyl group W{A 2 ) 
generating elements Wi,W 2 on the fundamental roots ai, 02 has the form 


WiCXi — —(Ti, W- 1 CX 2 — Oil “1“ 02') W 2 O 1 — Oi 02) W 2 O 2 — — 02‘ 

Using the formula ([5]) and the Cartan matrix 


= 


2 -1 

-1 2 


we rewrite the above relations for the fundamental weights 

tciAi = A2 — Ai, u'iA 2 = A2, W2\i = Ai, u' 2 A 2 = Ai — A2. 


( 22 ) 


(23) 


The action of the others group elements on the fundamental weights is determined by their 
representation by the generating ones 


ws = W1W2 , wa = W 2W1, W5 = W1W2W1, Wo = e. 


Moreover 

det wi = det W 2 = det = —1, 

and the determinants of the other elements of the Weyl group is equal to the unit. The Weyl 
vector has the form 

P = Ol + Q!2 = Ai + A 2 . 

Using these formulas and the following notations n = mAi + nX 2 , 4> = (po-i + '00^2; we 
represent the W (242)-invariant function of two variables <hm,n (E]) in the form 


^ _ 27rim0 27rin'i/> 1 27rim('i/t—</>) 27rin'0 1 27rim0 27rin(^—1/>) 1 27rim('i/t—</>) —27rin0 , 

I —27rim'i/> 27rin((j!)—-i/j) , —27rim'i/> —27rin(j!) 

I o o \ C' o • 


(24) 


In accordance with the formulas ® - (HU), we introduce the diagonal matrices 

All = diag(0, tjj - - (j), -'tp, -p)), 

A 2 = diag('0, -0, 0 - "0, -0, -0), 

and define the matrices Mi and M 2 

= A: = 1,2. 


( 25 ) 
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Then the function can be written in the following form 

Let us introduce the matrices 

R, = {I, - pM,)-\ R, = (Je - qM2)-\ 

where Jg is the unit 6 x 6-matrix and p, q are the real parameters. In view of the relations 


1 dR^ 


m\ dp 

the function <hm,n can be expressed as 


, 1 dRo 

= M™, — ^ 

p=o dq'^ 


= M", 


g=0 


*hm,n — 


1 d^+'^ti{RiR2) 


mini dp'^dq^ 
The matrix R 1 R 2 is diagonal and its trace has the form 

1 1 


p,g=0 


tr(i?i R 2 ) = 


(1 - (1 - pe2’"hb-'?i))(i _ 

1 1 
+ 4 


+ 


+ 


1 1 


(1 - pe-2’"'4’)(i - ge27ri(9i-p)) ('x _ ' 

Bringing the sum in the right hand side of this expression to a common denominator and ex¬ 
pressing the exponential coefficients before p, q by the complex-conjugate variables 2 ;, 4 defined 
by the formula 

^ z = ^$ 0,1 = (26) 

we obtain the generating function of the Chebyshev polynomials of the first kind associated 
with the Lie algebra A 2 


Fi„ = 


E Kijp^q^ 

i,j=0 


{1 — zp + zp'^ — p^)(l — zq + zq^ — q^) 
The coefficients Kij in (|27j) are given by the formulas 


(27) 


Koo 

= 6, 

Kio 

= -4^, 

K 20 

= 2z, 

Koi 

= -4z, 

Kn 

CO 

K 21 

= -2(T 

Ko2 

= 2^, 

Ki2 

= -2{z' 

K 22 

= zz — 


( 28 ) 
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The resulting function fl27p with coefficients fl28p coincides with the corresponding generating 
function for Chebyshev polynomials of the first kind, obtained in [20] (theorem 3.3). 

The polynomials calculated by fl27p are different in normalization from the ones introduced 
by Koornwinder in the works ra- If we put 

^0,0 ~ g *^0,0) — 2 ^m,0j ^0,n ~ 2 ^m,n ~ ^m,ny TTin ^ 0 , 

we obtain exactly the Chebyshev polynomials of the first kind Tm,n in the Koornwinder nor¬ 
malization. Let us list a few polynomials, taking into account that Tm,n = ^n,m; 


^ 0,0 = 1) 

^ 1,0 = z, 

^2,0 = — 2z, 

Til = zz — 3, 

T^fl = — 3zz 3, 

T 2 ,i = z‘^z-2z^-z), 

Tifl = z‘^ — Az^z + 2z^ — 2-1-42:, 
Ta^i = z^z — 3zz‘^ — z^ + bz, 

T 2,2 = Z^z^ - 2z^ - 2z^ + Azz - 3. 


The polynomials ‘hm,n satisfy recurrence relations of the form ([7]) 


^k+l,m, *^1,0 

^k,m+l *ho,l ^k—l,m ^k+l,m—l- 


(29) 

(30) 


These recurrence relation fl3UI) can be reduced to the following "linear" form by the simple 
transformation 


^k+l,m — *^* 1,0 ^k,m “ *^* 0,1 ^k-l,m + ^k-2,m, (31) 

‘hfc,m+l = "hop ^k,m-l + ^k,m-2- (32) 

Using the relations of a type similar to (l3Tp . fl3^ one can to obtain generating functions 
of Chebyshev polynomials by a scheme slightly changed from the present section. We will 
demonstrate it below in the section dedicated to the polynomials associated with the root 
system of the Lie algebra C 2 . 

2.2 The generating function for the polynomials of the second kind 

Now we calculate the generating function of the Chebyshev polynomials of the second kind 
associated with the root system of the Lie algebra A 2 . Taking into account the form of the 
Weyl vector p = Ai -|- A 2 we rewrite the general formula flT^ in the present case as 

wgW _^ 

X) dettne2^k«’P.</>) “ 

uisW 
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Un(4') 


(33) 




where the numerator and the denominator of (155]) are given by the 

following expressions 


— ('g27ri(m+l)<^g27ri(n+l)'!/’ _j_ g27ri(m+l)(i/)-0) g-27ri(n+l)0 _|_ g-27ri(m+l)i/)g27ri(n+l)(</.-p) ^ 
_ ^g27ri(m+l)(V>-0)g27ri(n+l)V> g27ri(m+l)</ig27ri(n+l)(</i-V>) _|_ g-27ri(m+l)V>g-27ri(n+l)<j!)^ 


= (^e27ri(</.+V>) g27ri(p-2<^) ^2ni{4>-2tp) ^ _ ^g27ri(2i/)-0) g27ri(20-p) g-27ri(0+V>) ^ ^ 

Let us introduce the following variables by 


Uw = 


<j)as 
^2,1 
<j)as ’ 
^1,1 


Uoi = 


$ 

¥ 


as 

1,2 


as 

1,1 


Using the relations fl55|) . (155]) . we obtain 


(34) 

(35) 


Ulo — z, Uqi — 


where and ^ are defined in (I55|) . 

In accordance with the method described in the Introduction, we define the following four 
diagonal matrices 

Mi+ = diag(e2"‘‘^, M 2 + = diag(e2-‘^, 

Mi_ = diag(e2"‘(^-'^), M^- = diag(e2"‘^, 

In terms of these matrices the function can be rewritten in the form 

= tr(M™ - M-M 2 U). 

Let us consider the matrices Ri± = {Iq —pMi±)~^ and R 2 ± = {lQ — qM 2 ±)~^. From the relation 


1 dR^^ 


m\ dp" 

it follows that the function 


p=0 


n\ dp"’ 


— M" 

— 2 ±) 


q=0 


K (p, q) = 


tr(i?i+i?2+ ~ R 1 -R 2 -] 


is the generating function of the two variables Chebyshev polynomials of the second kind 

1 d"^+"+^F^^{p,q) 


Um,.n 


(m + l)!(n + l)! dp"^+'^dq"+'^ 

After a calculation of the traces 

tr(i?i+i? 2 +) = ^ 

tr(i?i_i? 2 -) = 


p=0,q=0 


(1 — pe*^)(l — ge*-®) (1 — (1 — pe“®-®))(l — ’ 

111 

+ 


(1 — pe®(^“"^))(l — ge®^) (1 — pe®^)(l — ge®("^“^)) (1 — pe“®^))(l — ge“®A)' 
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and changing of the variables fl26p we obtain 


tr(i?i+i? 2 + ~ R1-R2-) 


pq{l— vq)^'i,i 

{1 — zp + zp^ — p^){l — zq + zq^ — q^) 


Dividing this expression on the singular element and on the product pq (to ensure that 
polynomial subscripts are in agreement with the order of differentiation) we obtain for the 
generating function of the Chebyshev polynomials of the second kind the known formula pO] 




Using 


(1 — Zp -\r Zp"^ 

we obtain by differentiation 


1 — pq 

p^)(l — zq + zq^ — q^) 


(36) 






m\n\ 


dp^dq^ 


p=0,(j=0 


(37) 


exactly the polynomials given by Koornwinder ini without additional normalization. We list 
a few polynomials, taking into account that Um,n = Un,m 


Uop = 1, 

Uip = 

f^2,0 = Z^ — Z, 

f/14 = zz — 1 , 

= z^ — 2 zz + 1, 
f/2,1 = z'^z — z"^ — z, 

= z^ — z 2z, 

Uz,i = z'^z - 2 zz^ - z^ + 2 z, 


3 The case of Lie algebra C 2 

3.1 The generating function for the polynomials of the first kind 

The root system of the Lie algebra C 2 has two fundamental roots ai, a 2 and includes the 
positive roots cxi + 02 , 2ai + 02 together with their reflections. Using the formula (jl]) and the 
Cartan matrix of the algebra C 2 


= ( -2 ~2 ) ■ 

we obtain the action of the Weyl group 1U(C2) generating elements wi,W 2 on the fundamental 
roots 

WiCni = —cti, Wicy .2 — 2 q,i T 02, 1^20^1 “ ^1 T 0^2, W 2 C (2 ~ —0^2- 
Next we rewrite these relations for the fundamental weights by means of the formula ([5]) 

rciAi = A2 — Ai, rciA 2 = A2, rc 2 Ai = Ai, W2X2 — 2 Ai — A2. 
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The action of the others group elements on the fundamental weights is determined by their 
representation by the generating elements 

W3 = W1W2, Wi=W2W3, W3 = W1W2W1, We = W2W1W2, WT = {wiW2f, 6 = Wq. (38) 
The determinants of the Weyl group elements are given by 


detwi = dettC2 = dettcs = dettcg = — 1, 
and the others are equal to the unity. The Weyl vector has the form 

3 

p = 2ai + -02 = Ai + A 2 . (39) 

Repeating the calculations from the Section 2.1 and using the notations which were intro¬ 
duced there we find W(C 2 ) - invariant function (orbit function) in two variables 


^ _ g 27 ri(m()i+np) _|_ ^ 27 ri(m(p-c/>)+np)) g 27 ri(m</i+n( 2 </)-p)) _|_ g 27 ri(m(V>- 0 )+n(- 20 +p)) 

_j_ g 27 ri(m( 0 -')/))-|-n( 2 </)-')/))) _|_ g 27 ri(-m(/)+n(- 2 </)+-!/))) g 27 ri(m( 0 —)/')-np) _|_ g 27 ri(-m 0 -np) 


(40) 


Note that the function iHOl) is real. The matrices Ai (|25|) in the considered case have the 
form 

Ai = diag(0, -0 - 0, 0, V' - 0,0 - V’, -0, 0 - "0, -0) 

and 

A 2 = diag(0,0, 20 - 0, -20 0, 20 - 0, -20 0, -0, -0). 

Then the function ^rn,n can be expressed by the matrices fc = 1, 2 in the trace 

form 

= tr(M0M2"). 

Following jlH], we define the Chebyshev polynomials of the first kind for algebra C 2 by the 
relations 


^ 0,0 — ^m,0 — RD,n — — ^m,n Tfl • fl ^ Q. 

O 2 2 

With this definition Tq^q = 1 and the new real variables x, y are given by the formulas 
a: = Ti^o = -t- 

y _ Tqi = + g27i-i(2?i-y) _|_ ^-2m{2<t>-p) ^ 


(41) 

(42) 

(43) 


Turning to the calculation of the generating function we introduce the matrices 

Rp = (-fs — pMi) Rg = {Is — qM2) 


which are diagonal with the elements given by 

(1 -pexp(27ri(Ai)i^J)"\ (1 - gexp (27^i(A2)kk))"^ 

where {Ai)^^ are the diagonal elements of the matrices Ai, i = 1,2, and /§ is the unit 8x8 
matrix. Substituting Rp, Rg in the equation (ITT)) and expressing the coefficients of p, q as 
functions of x (]T2|) and y (1T3|) . we obtain from (IT5l) 


= 

p,<i 


t,j=0 


(1 — xp -f (2 -|- y)p^ — xp^ -f p^)(l — yq + (x^ — 2y — 2)q^ — yq^ + q^) ’ 


(44) 


11 



where Kij are given by the formulas 

Kqo = 8, 

Kio = -6x, 

K20 = 4:1/ + 8, 

K30 = -2x, 

Km = -Qy, 

Kii = 5xy — 2x, 

K 21 = -4y‘^ + 2x^ - 10 ?/, 

Km = 2xy-2x, (45) 

Ko 2 = 4x'^ -8y -8, 

Ki 2 = —4x^ + 9xy + lOx, 

K22 = 8x^y — Qy^ + 4x‘^ — 20y — 8, 

K32 = —2x^ + 5xy + 6x, 

Ko3 = - 2 ?/, 

Ki 3 = 2xy - 2x, 

K23 = -2?/^ + 2x‘^ - Qy, 

K33 = xy - 2x. 

A few polynomials calculated using (144)) with the normalization (I4T)) are listed below 
^0,0 = 1 , 

^1,0 = X, 

T2,o = - 2?/ - 4, 

^3,0 = — 3xy — 3x, 

To,i = y, 

Ti,i = xy - 2x, 

T2 ,i = x'^y -2y^ -6y, 

Ts^i = x^y — 3xy^ — Axy + 2x, (46) 

To,2 = y^-2x^ + Ay + A, 

Ti ,2 = xy^ — 2x^ + 3xy + 6x, 

^ 2,2 = x^y^ — 2x^ — 2y^ — \2y^ + 8x^y + lOx^ — 29y — 8, 

^3,2 = x^y^ — 2x® — 3xy^ + lOx^?/ — 15?/^x + lOx^ — 2bxy — lOx, 

To,3 = y^ - 3x^?/ + 6?/^ + 9?/, 

Ti 3 = xy^ — 3x^y + 5?/^x + 2x^ + Qxy — 6x, 

T2,3 = x^y^ — 2y‘^ — 3x'^y — 16y^ + I2x^y‘^ + 29x^y — AOy^ — 30y, 

= x^y^ — 3x^y — 3xy^ + 15x^?/^ — 21xy^ + ISx^y — 45?/^x — 2x^ — 21xy + 6x. 

The recurrence relations for the polynomials under consideration can be obtained by the 
multiplication rules ([7]) putting s = (m, n) and k = (1, 0) in the hrst relation and k = (0,1) in 
the second one. As a result we hnally obtain 

‘hl,o‘hm,n = 4>m+l,n + 4*^-1, n + 4*^+1, n-1 + 4*^-1, n+1, (47) 

= ^m,n+l + ^m,n-l + ‘^’m+2,n-l + ‘^’m-2,n+l- (48) 
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For our purposes the relations (H7|) . fH5]) must be reduced to the "linear" form (see (1ST]) . 
(I32|) h Let us show the simple way to do this reduction. 

The determinant of the matrix pJg — Mi, gives us the characteristic equation, which satis- 
hes the matrix Mi. Since all the eigenvalues of this matrix have multiplicity 2, the minimal 
polynomial of Mi (after transition to the variables x, y) takes the form 

Pi = l — xp+ {2 + y)p^ — xp^ + p^. 

Therefore, the matrix Mi satisfy the equation 

Mf — xMi + (2 + y)M^ — xMi + /§ = 0. 

Multiplying this equation from the right by M("“^M;^ and taking the trace we obtain 

‘hm,n = X^rn-l,n ~ (2 + y)^m,-2,n + X^rn-3,n ~ (49) 

Acting in the same way with the matrix M 2 , we obtain the characteristic polynomial 

P 2 = l-yq+{x^ -2y- 2)q^ - yq^ + q^ 
and the second recurrent relation 


*hm,n — y^m,n—l (^ 2y 2)$m,n—2 T y^m,n—3 (^0) 

The linear relations (lT9l) and (l50l) can be tested independently by direct substitution of (H71) 
and (ITS]) . 

The relations (IT9|) and fl50|) allow us to hnd the generating function (ITT|) somewhat differently 
than it was done above. Let us rewrite them in the following form 


*hm,n ( ^m—l,n j ^m—2,n j ^m—3,n j ^m—4,n ) 

/ 
V 


( ^m,n—l ; ^m,n—2 i ^m,n—3 i ^m,n—4 ) 


/ X \ 

y + 2 

X 

V 1 

y 

x^-2y-2 

y 

1 


(51) 




(52) 


/ 


By means of the iteration of the relations (l5T]) and fl52|) it is not difficult to transform them to 
the form 


*hm,Ti ( ^m—l—k,n i ^m—2—k,n j ^m—3—k,n j ^m—4—k,n ) 

/ 

V 


*hm,n ( ^m,n—l—k j ^m,n—2—k j ^m,n—3—k j ^m,n—4—k ) 


f X \ 

y + 2 

X 

\ 1 
y 

x^-2y-2 

y 

1 


(53) 




(54) 
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where 


/ 


M. = 


X 


0 \ 


( y 

-1 

0 

0 \ 

0 

-1 

1 ~ 

x^-2y-2 

y 

0 

0 

-1 

0 

0 

-1 

0 


\ 1 

0 

0 

0 / 

= n 

- 4 in (l5Tp . 

We find 






( 

X 

\ 


^2,n 

$1,,; 

y + 2 

X 




tu—A 


\ 1 / 

/ y \ 

x^-2y-2 


V 


y 

1 


* hm,n — ( , ^ m ,2 ; *^^,1 ; *^^,0 ) 

Because of the matrices M^, My are inverted (det(Mj;) = det(My) = 1) we have 


/ X \ 


( y \ 


( 0 ^ 

y + 2 

= 

x^ -2y-2 


0 

X 

y 

y 


0 

\ 1 


\ 1 y 


v-iy 


Therefore, the relations ([53]), fl5T|) can be converted to the final form 


*hm,n — 


$ 


3,n ) 




2,n ) 


^1, 


-^o,n) m: 


( 0 \ 
0 
0 

v-iy 


(55) 


*hm,n — ( j 


-$ 


m,2 j ) 


-$r„,0 ) M, 


/ 0 \ 
0 
0 

V-1/ 


(56) 


The obtained relations (155 p and (15BP allows us to represent any Chebyshev polynomial ^m,n 
with subscripts m, n > 3 by the polynomials with 0 < m, n < 3. Indeed, using the relation (I55p 
with the polynomials ^rn,n, 0 < m,n < 3, one can to find the polynomials $m,o, ^m, 2 , 

‘hm ,3 with arbitrary subscript m. Then, using the relation (l56P with the calculated polynomials, 
one can to find a polynomial ^rn,n with an arbitrary subscript n. 

The considered modification of the method gives the advantage of the initial polynomials 
arbitrariness, that in turn allows us to construct different series of polynomials which satisfy 
linear recurrence relations (IT9l) and (l50p . but generally speaking not the original relations (IT7P 
and (jUj). In other words, the arbitrariness in the assignment of initial polynomials is limited. 
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3.2 The generating function for polynomials of second kind 

In the present Section we mainly follow the Introdnction and Section 2.2. Because of the Weyl 
vector for the algebra C 2 is equal to p = Ai + A 2 , the generating function for the Chebyshev 
polynomials of the second kind has the form (jHHj) 

TT (fk\ — _ 

n\(P) - 

uigW 

Using the notations of Section 2.2 and the presentation of the Weyl group IU(C' 2 ) elements by 
the generating ones of Section 3.1, we obtain for and the singular element respectively 


n+l 


= ^g27ri(m(/)+n'i/)) _|_ g27ri(m(p-(/i)+n(-2(/i+p)) _|_ g27ri(-m(/)-n'i/)) _|_ g27ri(m(0-'i/))+n(20-'i/i)) j _ 
_ ^g27ri(m(i/)-(^)+np)) g27ri(m(/)+n(2(/i-p)) g27ri(-m(/)+n(-2(/i+p)) _|_ g27ri(m((/i-p)-np)^ 

= (^g27ri(fli+p) _|_ g27ri(2p-3</)) _|_ g27ri(-(/)-p) g27ri(-2p+3fli) ^ _ 

_ ^g27ri(2p-0) g27ri(30-p) g-27ri(-3(/)+p) _|_ ^2m{<t)-2ij))y 

We not introduce at this stage new variables by the formulas 


(57) 

(58) 




45*1 

^1,1 


Uoi = 


<f,as 

^1,2 

(Jjas ’ 
^1,1 


because it will be more suitable to do later. Using the formulas fl57|) . fl58|) we define the following 
four diagonal matrices 

Mi+ = diag(e2^‘‘^, g2^i(<A-h))^ M2+ = diag(e2^‘^, g- 2 -ip^ g2-i(20-p)^ 

Mi_ = diag(e2^‘(’^-'^), M2- = diag(e2^‘^, g2^i(20-P)^ g2^i(-20+p)^ g-27rip^_ 

In terms of these matrices the function has the form 

= tr(M- 

Introducing the matrices Ri ± = (J4 — i?2± = (-^4 — ?d^2±)~^ one can to say that the 


function 


K (p, q) = 


tr(i?i_|_i?2+ ~ R 1 -R 2 -] 

*^ifi 


is the generating function of the two variables Chebyshev polynomials of the second kind 

1 d^+^+‘^F^^{p,q) 


(m + l)!(n + 1)! dp"^+'^dq^+'^ 


p=0,q=0 


As a result of the simple calculations we obtain 

pg(l + q + p^q + p^q^ — xpq)^\\ 


tr(i?i_|_i?2+~-R 1 --R 2 -) — 


{1 — xp+ {y + 2 )p 2 — xp'^ + P^)(l — yq + — 2|/ — 2)q‘^ — yq^ + q^) ’ 
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where was introduced the following intermediate notation 


X = tr(Mi+) = tr(Mi_) = 

y = tr(M 2 +) = tr(M 2 _) = + e2-i(20-V>. 

Now as in the A 2 case we divide the obtained expression by the singular element and by 
the product pq. The calculation of derivatives gives us 


Uio = 


dF;\p,q) 


dp 


= X, 


Uoi = 


dFy{p,q) 


p=0,q=0 


dq 


= y + l. 


p=0,q=0 


Just in this place it is convenient to introduce the new variables 


X = X, y = y + 1. 

In the terms of these variables F^\p, q) has the form 

\ + q + p^q + p^q^ — xpq 


F (p, q) = 


{1 — xp + {y + l)p^ — xp^ + p^)(l — [y — F)q + — 2y)q^ — {y — 1)?^ + (f") 

(59) 

Using the formula 

2 {p,q) 


Um..n 


min! dp'^dq^' 
we calculate a few polynomials of the second kind 


p=0,q=0 


Uo,o — 1, 

U\fl = X, 

Uo,i = y, 

U2,o = x'^-y- 1 , 
f/1,1 = xy -X, 

Uo ,2 = -x'^ + y‘^ + y, 

Uzfl = x^ - 2xy - X, 

U2,i = x^y - x^ - y‘^ - y+ 1, 

Ui^2 = —x^ + xy"^ + x, 

Uo,3 = - 2x^2/ + + 2y‘^ - 1. 

These polynomials coincide with the ones which were found in [11] from the recurrence relations 
if interchange the polynomial subscripts and perform the following substitution x ^ X 2 , p —)■ 
Xi. 


4 The case of the Lie algebra G 2 

In this Section we calculate the generating function of the Chebyshev polynomials of the first 
kind associated with the exclusive Lie algebra G 2 - 
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The set of the positive roots of the algebra G 2 root system includes two fundamental roots 
«!, a 2 together with ai + 02 , 2ai + Q! 2 , 3ai + Q! 2 , 3ai + 2a2. Using the relation (jl]) and the 
Cartan matrix 


Cg, 


2 

-3 2 r 


we obtain the action of the Weyl group W(G 2 ) generating elements Wi,W 2 on the fundamental 
roots 

W\(y.\ = —tti, WiO'2 ~ 3q;i T 0'2i ^ 20^1 ~ T tt2, W2C(2 ~ — 0 ^ 2 - 


Taking into account (|5]) we obtain for the fundamental weights 


rciAi — A2 — Ai, W1X2 — A2, W2X1 — Ai, W2X2 — 2 Ai — A 2 . 


The simple calculation results in the following W-invariant functions corresponding to the group 

W(G2) 

(|) _ g27ri((2m+n)<^+{m+2/3n)i/>) _j_ g— 27ri((2m+n)0+(m+2/3n)p) _|_ g27ri((m+n)</)+(m+2/3n)i/>) _j_ 

g-27ri((m+n)<ji+(m+2/3n)'i/)) _|_ g27ri((2m+n)(ji+(m+l/3n)'i/j) _|_ g—27ri((2m+n)</)+(m+l/3n)-(/') _j_ g27ri((m+n)(/)+l/3n'i/') _j_ 
g—27ri((m+n)<^+l/3np) _j_ g27ri(m<^+(m+l/3n)'i/j) _|_ g—27ri(m0+(m+l/3n)-(/)) _j_ g27ri(m0—l/3n'i/)) _|_ g—27ri(mi/)—l/3np) 

(60) 


This formula is consistent with the orbit function which was obtained in |18) . 

Let us introduce the diagonal matrices 

Ai = diag(20 + V’, -20 - 0, 0 + 0, -0 - 0, 20 + 0, -20 - 0, 0, -0, 0 + 0, -0 - 0, 0, -0), 

A 2 = diag(0 + 2/30, -0-2/30, 0 + 2/30, -0-2/30, 0 + 1/30, -0-1/30, 0 + 1/30, -0-1/30, 

+ 1/30,-1/30,1/30,-1/30). 

Then, as in the case of the algebra G 2 , we dehne 

Mfc = fc = 1, 2, = tr(M0M2”). 

We introduce by dehnition the new real variables 


a: = tr(Mi)/2 = 


_|_ g-27ri(20+i/)) g27ri(0+i/)) g-27ri(0+i/>) _|_ g27ri0 g-27ri(/.^ 


(61) 


and 

y = tr(M2)/2 = + e“^’^‘U+2/3V') _j_ g27ri(0+l/3V’) _|_ g-27ri(0+l/3V’) _j_ g2/37riV> _|_ g-2/37riV'_ 

Since all eigenvalues of the matrices Mj, i = 1,2 in the considered case have multiplicity 2, 
then the degree of its minimal polynomial is equal to 6 

Pi = 1 — xp + {y^ — 3xy — 5x — 9y — 9)p^ — {—2y^ + + 6xy + 12x + 18y + 20)p^ + 

+ {y^ — 3xy — 5x — 9y — 9)p^ — xp^ + p®, 

P 2 = I — yq + {y + X + 3)q^ — {y'^ — 2x — 4)q^ + {y + x + 3)g"^ — yq^ + g®. (63) 
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Using the method described in Section 3 we obtain a linear recurrent relations 


<hm,n - + {y^ - ^xy -5x-9y - 9 )^m-2,n- 

{-2 y^ + x^ + 6 xy + 12a; + 18y + 20 )<l)m- 3 ,n+ 

{y^ - 3xy -5x-9y- 9)^rn-4,n - X^m- 5 ,n + ^m-6,n = 0, 

*hm,n y^m,n—l~\~{yj~\~X-\-3')^rn,n—2 2 x (?/-t-X-|-3)$m,,n—4 6 ~ 0- 

Acting by the above technique, we obtain the desired generating function in the form 

F,., = (PiA)-' ( ^ A«pV ) . (64) 

\i,i=o / 

where the polynomials Pi, P 2 are given by the relations (l62|) . ([63]), and the coefficients Kij have 
the form 

A'oo = 12, 

Pio = -lOx, 

-P 20 = = 8?/^ — 2Axy — 40x — 72y — 72, 

A’so = —6x^ + 12y^ — 36xy — 72x — 108?/ — 120, 

-P 40 = 48?/^ — 128x8?/ — 208x — 368?/ — 36, 

K 50 = -2x, 

Kqi = - 10 ?/, 

Pii = 9xy + 2y — 2y‘^ + 4x + 12, 

K 21 = 40x?/ — 8?/*^ + 25x?/^ — 2x^ + 72?/^ — 6x + 74?/, 

K^i = Qx^y — 12y'^ + 35x?/^ + 2x^ + 72xy + 108?/^ + 6x + 118?/, 

K 41 = —4?/'^ + 12x?/^ + 21xy + 38?/^ — 4x + 34?/ — 12, 

Psi = 2 xy - 2y, 

Ko 2 = 8x + 8?/ + 24, 

Ki 2 = —I3xy — 8x^ — 16?/ + 2y^ — 28x — 12, 

K 22 = Qy^ + 7xy^ — 21x^y — 18x?/^ — 40x^ — 166x?/ + 24?/^ — 56?/^ — 192x — 268?/ — 216, 

K ^2 = 13x?/^ — 44x^?/ — 6x^ +10?/^ — 92x^ — 29xy‘^ + 36?/^ — 285xy — 90?/^—422?/ — 338x — 348, 

-P 42 = —I2x^y + 4x?/^ — 13xy^ + 4?/"^ — 92x — 20x^ — 128?/ — 86 xy — 36y^ + 10?/^ — 96, 

K 52 = —2xy — 2x^ + 2 y‘^ — 2y — 6x, 

Kqs = —6?/^ + 12x + 24, 

Ki 3 = 6xy^ — 2xy — 12x^ — Qy — 26x, 

K 23 = —Qy^ — 3bx^y + 39xy^ + 31xy‘^ — 62x^ — 187xy + 82?/^ + 54?/^ — 254x — 246?/ — 252, 

7733 = 5x^?/^ — 19y^ — lOx^ + 60x?/^ + 50x?/^ — 62x^?/ — 339xy + 138?/^ — 144x^ + 

+ 96?/^ - 432?/- 480x - 456, 

7743 = —4?/^ + 20x?/^ — 23x^?/ + 21xy‘^ — I27xy + 56y^ — 42x^ + 36?/^ — 174?/ — 178x — 180, 
7753 = 2 xy‘^ — 2xy — 4x^ — Qy — lOx, 

77o 4 = 4x + 4?/ + 12, 
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i^i4 = —Axy — Ax^ + 2y‘^ — 2y — 12x, 

K2A = —I2x^y + Axy^ — 13a:?/^ + Ay‘^ — 92x — 20x‘^ — I28y — 86xy — 36y^ + 10?/^ — 96, 

= —SOx^y—4x^ — 22Qx + Gy"^ — 1 7xy'^ + 9xy^ — 189xy — 62x^ + 2Ay^ — 5%^ — 2,7 Ay — 228, 
Kii = 3xy^ — 9x^y + 2y'^ — 6xy^ — 7Axy — 20x^ + 12y^ — 20y^ — 96x — 12Ay — 108, 

-ft '54 = —2x‘^ + 2y^ — 7xy — lOo: — 16y — 12, 

Ko 5 = —2y, 

Ki 5 = 2xy - 2y, 

K25 = —2y‘^ + Gxy'^ + llxy + 29y‘^ — Ax + IGy — 12, 

/^35 = 2x^y — Ay* + llxy"^ + 2x^ + 2Axy + 36?/^ + 6x + 38y, 

-ft'45 = —2y* + 7xy‘^ — 2x^ + 19xy + 18y‘^ — Gx + 20?/, 

-^55 = xy + 2y — 2y‘^ + 4x + 12. 

In conclusion we list here a few polynomials generating by the function flM|) 

^00 = I5 ^10 = X, Toi = ?/, 

T 20 = — 2y^ + Gxy + 18y + lOx + 18, 

Til = xy — 2y‘^ + 2y + Ax + 12, 

To 2 = y"^ -2x-2y-G 

T 30 = x^ — 3y^x + 9x^y — Gy^ + 18a;^ + Abxy + Gy + G3x + 60 
T 21 = x^y — 2 y* + Gxy"^ + 2 x^ + 12 xy + 18y‘^ + Gx + 20 y 
Ti 2 = xy'^ — 2x^ + 2y‘^ — 3xy — lOx — Ay — 12 
To 3 = - Gxy -Gx-9y- 12. 

By the method described in the previous Sections we can obtain the generating function of 
the Chebyshev polynomials of the second kind, associated with the algebra G 2 - Here are a few 
polynomials calculated in the work [19] from the recurrence relations 

Uqo = 1 , UiQ=x, Toi = ?/, 

U20 = 2y + x A-x"^ - y^ + 2xy, 

Uii = 1-y^ + X + xy, 

U 02 = -l-x-y + y'^, 

U30 = —y'^ + 2x + 3x^ — y^ — 2xy^ + y* + x^ + Axy — 2xy‘^ + Ax^y, 

U21 = -l-y + 2 y‘^ + x^ + y^-y* + xy‘^ + x^y, 

U 12 = y+ y‘^ - X - x'^ -y^ + xy"^, 

Uo3 = -y-y^ -x + y^ - 2 xy. 
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